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I. INTRODUCTION 

In his seminal paper [1], Werner Israel considered the 
dynamics of a self-gravitating thin matter shell. The 
main purpose of his theory was to find a simple model 
describing gravitational collapse. In case of a realistic col- 
lapse, equations describing evolution of matter and grav- 
ity are extremely difficult to handle. Israel's idea was 
that many aspects of the collapse may be investigated 
within a toy-model, which consists of a matter shell and 
the surrounding gravitational field. The dynamics of such 
a system reduces to a proper tailoring of the two different 
vacuum solutions describing the two sides of the shell. 

We are going to present a systematic derivation of the 
Israel model from the variational principle and the con- 
struction of its canonical (Hamiltonian) structure. In 
our approach the space-time M consists of two parts, 
tailored together along a hypersurface S, which contains 
a moving matter shell. "Tailoring" means that the in- 
duced metric g a b of S is continuous. On the other hand, 
the four-dimensional connection coefficients may be 
discontinuous on S. It is proved that the singular part 
of the Einstein curvature tensor density of such a space- 
time contains derivatives of those discontinuities and may 
be defined in the sense of distributions as G a b — G a i,ds, 
where 5s is a Dirac delta distribution concentrated on S. 
As we show in the sequel, the following relation holds: 

G\ = [Q\] , (1.1) 

where square brackets denote jump of the extrinsic cur- 
vature across S, written in the ADM form. The sin- 
gular Einstein tensor density matches the singular (con- 
centrated on S) energy-momentum tensor density of the 
matter shell. Distributional Gauss-Codazzi equations 
are then derived (and not postulated, as in the Israel 
approach). They imply that both the Einstein tensor 
density and the matter energy-momentum tensor den- 
sity are conserved. We prove that the conservation 
law V a G a b = can be written in terms of the three- 
dimensional geometry of S. 



The possibility of defining the singular Einstein tensor 
and calculating its divergence (in the sense of distribu- 
tions) via the standard formulae of Riemannian geome- 
try, simplifies dramatically the calculational part of the 
theory. For this purpose we must assume that the four- 
dimensional metric is continuous across S. At this point 
our approach differs from the technics used by many au- 
thors, including W. Israel himself, who always stressed 
that no assumptions about the continuity of the four- 
metric, except for the continuity of the three-metric on 
the surface S, are necessary. This is, of course, true. We 
observe, however, that our "additional condition" on the 
continuity of the entire four-dimensional metric, does not 
contain any geometrical or physical condition imposed 
on configurations of the physical fields considered by us, 
but is merely a gauge condition imposed on coordinate 
systems used. Whenever an intrinsic three-dimensional 
metric of S is continuous, then also the remaining four 
components of the four-dimensional metric may become 
continuous after a suitable coordinate transformation. In 
this new coordinate system we can use our technics based 
on the distribution theory, but the dynamics derived this 
way is written in terms of intrinsic, geometric relations, 
having sense in an arbitrary system of reference. Hence, 
our derivation does not depend upon particular gauge 
conditions which we use. Note that even in a perfectly 
smooth, flat space-time, we can introduce a coordinate 
system, in which only intrinsic three-dimensional metric 
on a fixed hypersurface S = {x 3 = const.} is contin- 
uous, whereas the other four components of the four- 
dimensional metric g^^ may be discontinuous. Using 
these "singular coordinates" one can properly formulate 
e.g., the initial value problem for the Maxwell field, but 
nobody uses such a formulation for obvious reasons. Our 
additional assumption on the continuity of space-time 
metric is motivated by similar reasons: not to compli- 
cate things which already are complicated enough. The 
technical flexibility we obtain this way, i.e., the possibil- 
ity of defining the singular geometric objects living on S 
via standard geometric formulae, where the derivatives 
are understood in the sense of distribution, is the main 
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profit of our approach. 

The structure of the paper is following. At the be- 
ginning we derive Israel's model from variational princi- 
ple, first in the Lagrangian formulation, and then - after 
performing an appropriate Legendre transformation - in 
the Hamiltonian picture. Many technical elements of the 
construction presented here come from the papers [2] and 
[3]. A significant simplification of the calculations is ob- 
tained using consequently tensor densities instead of ten- 
sors for such quantities like the matter energy-momentum 
and the Einstein curvature object (see also [4]). Part of 
presented calculations was included in PhD thesis of one 
of us [5]. 



II. PROPOSED DESCRIPTION OF THE 
ISRAEL'S MODEL 

Consider a space-time consisting of two parts which 
are tailored together along a hypersurface S, whose non- 
degenerate metric has the signature (—,+,+). Unlike in 
the original Israel's approach (which was also used in [2]), 
here we restrict ourselves to coordinate systems for which 
all the components of the space-time metric are contin- 
uous. As was already mentioned, this condition does 
not limit the applicability of our formalism and may be 
treated as merely a " gauge condition" imposed on the co- 
ordinate system. It simplifies dramatically theoretical de- 
scription of the model. Further simplification is obtained 
by using a coordinate system such that the hypersurface 
S is given by the equation {x 3 = const.}. Metric deriva- 
tives along S (i.e., d a g^v, where a — 1, 2, 3 ) are continu- 
ous, whereas the transversal derivative (i.e., d^g^v) may 
have jumps over S. We assume that the topology of £ 
is of the type S 2 x R 1 , i.e., it describes a history of a 
matter concentrated on a two-dimensional surface with 
the topology of the sphere S 2 . 

The shell divides the space-time into the internal and 
the external part with respect to the world tube S. In 
both parts vacuum Einstein equations may be derived in 
a standard way from the variational principle. Hence, 
the regular part of the Einstein tensor must vanish ev- 
erywhere outside and inside of S. Only its singular part 
concentrated on S is left. 

The singular part of the Riemann tensor is propor- 
tional to the (invariant) Dirac delta distribution 5s con- 
centrated on S, because first derivatives of the metric 
(and, whence, also the connection coefficients ) may 
be discontinuous across S. In our particular coordinate 
system we have 5 s = 5(x 3 ). As will bee seen in Section 
IV, the resulting singular part of the Einstein tensor will 
match the singular energy-momentum tensor describing 
matter concentrated on S. This singular part may be 
obtained from the standard formula for the Ricci tensor: 



following combination of Christoffel symbols: 
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Then we have: 
R^iv — 



A x A a 



— A x A a 

[lX VG ' 



(2.2) 



(2.3) 



Because A may have only discontinuities of the "jump- 
type" across S, the derivatives of A along directions tan- 
gent to S are thus finite and belong to the regular part 
of the Ricci tensor. Hence, its singular part consist only 
of the transversal derivatives: 



smg(R^) = d 3 A* = 5(x 3 )[A 3 } 



(2.4) 



where the square brackets denote the jump of a specific 
quantity across S. Hence, we have the following formula 
for the singular part of the Einstein tensor density: 
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sing := y/\g\ sing {R^--R\= 5(x 6 )G^ 



where 



ft 



[Alp] 



(2.5) 
(2.6) 



is a quantity living on S. Now, we are going to show that 
it is actually a 3-dimensional tensor density on S. For this 
purpose we first show that its components transversal to 
S vanish and, whence, its only non-trivial components 
are those tangential to S. To prove that it transforms 
like a tensor density on S let us observe that Q^ v be- 
haves like a 4-dimensional density, which splits into the 
3-dimensional density on S and the one-dimensional den- 
sity along x 3 . But the Dirac delta S(x 3 ) is already the 
density (and not a scalar!) on the real axis x 3 which 
proves that the remaining object G behaves indeed like 
a 3-density. Hence, we have to prove the following 



Lemma 1. 



G\ = 0. 



(2.7) 



Proof. On both sides of S consider the following combi- 
nation of the connection coefficients: 



Q% ■= M (gTAl v \^A 



(2.8) 



It is useful to encode the entire information about the 
metric in the following tensor density 1 
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(2.9) 



R^v — dxT* v - 9( M r^ )A + r* A r^ - r^ CT r^ A . (2.1) 

A considerable simplification is obtained if we use the 



1 In the so called affine formulation of general relativity, proposed 
by one of us in 1978 (cite [6]), this quantity plays role of the 
momentum canonically conjugate to connection T. 
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Hence, we have 

^Q"„ = ^ a A\ v \5^A% , (2.10) 
Then the equation (2.6) gives 

G"„ := [Q"„] . (2.11) 

Now, we use metricity condition for the connection T, 
which is fulfilled on both sides of S (indices a, b = 0, 1, 2 
label the coordinates on S): 

= V a ^ 33 = d a ir 33 + 2-n 3 »T% - n 33 T^ 

= cU 33 + 27r 3 <M 3 a (2.12) 

o = v a ^ 3a = ^ 3a + 7r^r 3 a + ^r; a -^ 3a r^ 

= d a 7r 3a +7r ab Al b -7r 33 A 3 3 , (2.13) 

But the metric components 7r Miy and their derivatives 
along S are continuous across S. Hence, we obtain 



±-G 3 a = ^[Al a ] = -\[d a , 33 ]=Q 

^G 8 , = ~\ {? ah \A\ h \ - , 33 [A%,\) = \[d a ^] = , 

(2.15) 



which ends the proof because in our coordinate system 
we have G^ v = G 3 „. □ 



III. GEOMETRIC INTERPRETATION OF THE 
SINGULAR CURVATURE AND GENERALIZED 
BIANCHI IDENTITIES 

Equation (2.11) expressing the singular part of Ein- 
stein tensor in terms of jumps of quantities Q^ v across 
S is extremely useful in our derivation of the shell dy- 
namics. However, is not satisfactory from the geometric 
viewpoint because quantities Q^ u do not have any ge- 
ometrical meaning on both sides of S, and only their 
jump (2.11) across S does. Now, we are going to prove 
that we obtain the same result replacing non-geometric 
object Q^ v by a tensor density Q^ v which is, by defini- 
tion, orthogonal to S and whose restriction Q a b to S is 
equal to the external curvature of the hipcrsurface S in 
the ADM ([7]) representation. In our coordinate system 
S = {x 3 = const.}, external curvature of S is given by: 



L a b ■' 



■ ab 



l ab i 



(3.1) 



and its ADM version equals (cf. [8], [3]): 

Q ab := Vldct^l {Lg ab - L ab ) . (3.2) 

Here, g ab is a three-dimensional inverse of the induced 
metric g a b on S and L ab = g ac L c dg db - It is easy to 



check that g ab may be calculated in terms of the four- 
dimensional inverse metric g^ u via the following formula 
(see [8]): 



„3a 36 

n ab = n ab - y y 

y y n 33 



(3.3) 



In order to express coefficients Q^ v in terms of the tensor 
density Q^ v , observe that identities (2.12) and (2.13) can 
be solved algebraically with respect to A\ z and A\ a . As 
a result, we have on both sides of the hipersurface S the 
following identities : 

Al^^(d a ^ a + Al b ^ b ) , (3.4) 
A 3 a = -^(d a 7r 33 + 2A 3 ab n 3b ) . (3.5) 

Hence, all the coefficients A 3 ^ can be expressed in terms 
of A 3 b , i.e.,, using (3.1) and (3.2), in terms of Q ab and 
the metric. Using again formula (3.3) we obtain: 



(2.14) Q» b = Q a b+ 167T 



1 7T 3a 1 / 7T 3C 

1 71 ^33 , 1 ra ( 71 _33 _3c 
7T h + 7:0 h 



2tt 33 ' b 2 b \ 7T 33 < c 



7T „ — 7T 



(3.6) 

The last term is identical on both sides of S, because 
the metric components are continuous. Hence, their 
jump across S vanishes and, due to (2.11), we obtain: 



G a b := [Q\ 



(3.7) 



whereas the transversal components G\ vanish. Because 
the object Q a b is a well defined tensor density on both 
sides of S, (its definition does not depend upon coordi- 
nate system used), hence the tensorial character of the 
three-dimensional object G a h has been proved. 

Now, we are going to show that the total Einstein ten- 
sor Q^ v of our spacetime M fulfills Bianchi identities. Be- 
cause regular part of Q^ v is discontinuous across S and, 
moreover, it contains also a Dirac-delta-like singular part, 
these identities must be understood in a distributional 
sense. To prove them we shall use Gauss-Codazzi equa- 
tions, relating transversal components Q\ of the Einstein 
tensor density Q% to a divergence of external curvature 
Q on S: 



< b 



V a Q a b = . 



(3.8) 



where V denotes the intrinsic, three-dimensional covari- 
ant derivative on S. But the transversal component Q\ 
is a well defined three-dimensional object on S. In our 
coordinate system, adapted to S in such a way that the 
coordinate a; 3 is constant on S, this quantity is simply 
equal to a "third" component: Q\ = Q\. Taking a jump 
of this equation across S we obtain the following identity: 



MQ)\\ + v a G a b = o 



(3.9) 



Now, we are going to show that the left hand side 
is exactly the singular part of the Bianchi divergence: 
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\7 = 0. Indeed, the regular part of this quantity 
vanishes on both sides of S (a consequence of the stan- 
dard Bianchi identities), whereas its singular part is pro- 
portional to 5s- To prove this statement observe that 
Q\ = and, whence, no derivative of the Dirac delta 
is produced when we apply the covariant derivative V3 
to the singular tensor density (2.5). What remains are 
thus the "along 5" derivatives V a . As a result of this 
operation we obtain, therefore, the quantity V G a 6 mul- 
tiplied by 8s- Another £-like term is obtained from the 
regular part reg(<?) which is discontinuous across S. Tak- 
ing the derivative V3reg(C7) 3 b of the regular part we ob- 
tain, therefore its jump of it across S multiplied by the 
Dirac delta: [rcg(Q) 3 b ]Ss- Finally, the singular part of 
the Bianchi identities is the sum of the above two ex- 
pressions: 

V M 0" 6 - ([reg(S)i] + V a G a b ) d(x 3 ) = , (3.10) 

where the last identity is just the Gauss-Codazzi equation 
(3.9). Hence, we have shown that identities V M <7 M fe = 
are also fulfilled for space-times with a singular curvature. 



metric tensor g^ and to the matter fields z K living on S. 
At the moment we do not specify the nature of the matter 
fields. It is enough to assume that they are geometric 
objects living on this surface and that matter Lagrangian 
L mat is a scalar density on S, which depends locally on 
the values of those fields, their derivatives along S and 
the metric of the surface S. 

From the point of view of the two regular half- 
spacetimes (which are tailored across S) the singular part 
of the action arises as the sum of the boundary contribu- 
tions from both the sides of the shell. On the other hand, 
the action (4.3) does not contain any surface term at in- 
finity. This is due to the techniques used here (cf. [3]), 
where we first derive the field dynamics within a spatially 
compact region D and then shift its boundary dD to the 
space infinity. Of course, the boundary manipulations at 
infinity are still necessary but in our approach they arise 
as a Legendre transformation between different control 
modes at the boundary dD (see Section V A) . 

There are many ways to calculate variation of the 
Hilbert Lagrangian. Here, we use a method proposed by 
one of us (see [10]). It is based on the following, simple 
observation: 



IV. DYNAMICS OF THE "SHELL + GRAVITY" 
SYSTEM 

Dynamical equations of the physical system composed 
of a matter shell and the surrounding gravitational field 
will be derived from the action principle 5 A = 0, where 



1 

16tt 



(4.4) 



A = AZ 



ATI + A B 



(4.1) 



is the sum of the gravitational action and the matter ac- 
tion. Gravitational action, defined as the integral of the 
Hilbert Lagrangian L = j^^/gR, splits into the regular 
A r ° r l v and the singular part A"^ , according to the decom- 
position of the curvature R = reg(i?) +sing(i?) (a similar 
"mixture" of a "bulk action" and a singular "body ac- 
tion" concentrated on a submanifold was recently used 
by C. Barrabes and W. Israel - see [9] - to derive brane 
dynamics in general relativity). 

Using formulae (2.5)-(2.7), we express the singular 
part of R in terms of the singular part of the Einstein 
tensor: 



Expressing R^ v in (4.4) by the connection coefficients 
L* and their derivatives, it is easy to show that the last 
term on the right hand side is a complete divergence due 
to the following identity: 



w^sr^ = d K (*r K Kv) = (W K ) K» 



(4.5) 



where, besides of the quantity (2.9), we have introduced 
the following notation: 



x , 



(4.6) 



and 



mnL-Z = V\9\ sing(i?) = -sing(S) 
= -G^ g ^S(x 3 ) = -G ab g ab 6(x 3 ) . 

Hence the total action is the sum of three integrals: 



(4.2) 



A 



D 



L rog 

grav 



D 



Dns 



(4.3) 



where D is a spatially compact four-dimensional region 
with boundary in M, which is possibly cut by a three- 
dimensional surface S (actually, because of the Dirac- 
delta factor, the second term reduces to integration over 
DC\ S). Variation is taken with respect to the spacetime 



d K T 



(4.7) 



Proof of the identity (4.5) is straightforward if one uses 
the fact that are not independent quantities but the 
connection coefficients, i.e.,, combinations of the metric 
components g^ v and their derivatives. This means that 
the covariant derivative \7ti with respect to connection T 
vanishes identically. Taking into account the fact that tt 
is a tensor density we have that: 



K \ — 1 Xk ^X 1 an 



(4.8) 



which implies identity (4.5) in a simple way see [10]). 
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Finally, variation of the regular part of the Hilbert La- 
grangian is following: 



A. Matter Lagrangian and the singular 
energy- momentum 



1 



\g\ reg(i?) ) = -— veg(grdg^ 
f reg(d K (n ) r K <*)) ■ 



\i67r v ,c " ' ° x ~' j i67r " N ' (4.9) 



We are going to show now that the same equation holds 
for the singular part of R, i.e., 



1 



<*0 = M^Vlsl sing(i?) 



16tt 



sing^)^^ + sing(a K (^ K ^)) . 

(4.10) 



Proof of (4-10): Calculate the singular part of 
® K {' 1T \ tVK ^^iv) ■ Because all these quantities are 
invariant, geometric objects (5T is a tensor!), we may 
calculate them in an arbitrary coordinate system. 
Hence, we may use our adapted coordinate system, 
where coordinate x 3 is constant on S. Taking into 
account the continuity of -k^ uk across S we obtain 

sing (d K (^KJ) = S(x a )nr ± 5[T^} 
= S(x^r 3 S[T^]=S(x^S[AlJ. 

From the definition (2.10) of the object Q^ v it follows 
immediately that 



J SA 



Using (2.11) we have that 

[Q^]S(x 3 ) - sing(g)^ . 
Putting these formulae together we obtain 



(4.12) 



(4.13) 



1 



6(x*)nr 5FU = singer 

= SL T™ + ^m g (grs 9fl „ , 

which ends the proof of (4.10). 



(4.14) 



□ 



Summing up the regular part (4.9) and the singular 
part (4.10) we obtain variation of the whole gravitational 
Lagrangian: 

SL s ^ = -^-g^S 9liV + d K (-K x ^ K ST^) , (4.15) 



Now, consider the matter part of the action. As was al- 
ready mentioned, we assume that all the dynamical prop- 
erties of matter fields are described by a Lagrangian den- 
sity L me , t , which is an invariant, three-dimensional scalar 
density on S. It depends upon some matter fields z K 
living on S, their first derivatives along S: z K a := d a z k , 
and the metric tensor g a b on S: 



L mat (z K ; z K a ; g ab ) ■ 



(4.16) 



Calculate the variation of the matter Lagrangian: 



SL^ t = 



<9L„ 



dg a i 



'-Sgab + 



dL r 



^■Sz K 



dz K 

2 T &9ab+ \~d^ 



dL mat 
' dz K a 

dz K a 



d a Sz 

5z K 



K 



-8a ( P K a Sz 



(4.17) 



where we have introduced the following three- 
dimensional symmetric energy-momentum tensor 
density T ab on S: 



T ab ._ 2 



dL 
dg ab 



(4.18) 



and the momenta px a canonically conjugate to material 
variables z K : 



PK ■= 



dz* a ■ 



(4.19) 



Finally, we obtain the following formula for the variation 
of the total ( "matter + gravity" ) Lagrangian: 

5L = - ±-r G g(gr5g uv 

, a. ( <9L mat dL 

mat \ ? K 



-6(x 3 )— (G ab ~8TTT ab ) Sg ab 

+ 8 K (*r K K») + S ( x3 ) d " (PK a 5z K ) . (4.20) 

In this Section we assume that both Sg^ and 8z K 
vanish in a neighborhood of the boundary dD of the 
space-time region D (this assumption will be later re- 
laxed, when deriving Hamiltonian structure of the the- 
ory). Hence, the last two boundary terms of the above 
formula vanish when integrated over D. Vanishing of 
the variation SA — with fixed boundary data implies, 
therefore, the Eulcr-Lagrange equations for the matter 
field z K , together with Einstein equations for gravita- 
tional field. Regular part of Einstein equations: 



which generalizes the corresponding formula from the pa- 
per [3] to the case of spacetimes with a singular curvature. 



reg^f = 



(4.21) 
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is satisfied in the whole space-time, whereas the singular 
part 



G ab = ^ab 



(4.22) 



must be fulfilled on S. This way the Israel equations for 
the shell dynamics have been derived from the variational 
principle. 

Information encoded in the three-dimensional object 
r ab may be also expressed in terms of a four-dimensional 
energy-momentum tensor T' 1 " := 8{x )t^ v , where, in 
agreement with (2.7), the transversal components of r 
vanish by definition: t ± " = (in our adopted coordinate 
system it simply means that t 3i/ = 0, but this condition 
in its previous form does not depend on the coordinate 
system). Summing up the singular and the regular parts 
of the gravitational field, we write the variation of the 
total Lagrangian as 

1 



16tt 

/ dL 
+ d(x 3 ) 1 



d a 



dL„ 



5z K 



\ dz K a dz K 



(4.23) 



which is the starting point of our derivation of the dy- 
namics of the system. We stress that this equation is an 
identity, implied by the structure of the action (4.1). 

Field equations of the theory, i.e., Einstein equations 
for gravitational field and the Euler-Lagrange equations 
for matter fields, are equivalent to vanishing of the vol- 
ume part of the above variation (i.e., of the first two 
terms): 



dz K 



8nT^ , 
p. <9L mat 
a dz^ a - 



(4.24) 
(4.25) 



This, in turn, is equivalent to the fact that for arbitrary 
variations (non necessarily vanishing on the boundary!) 
of the independent fields Sg^ and 5z K , variation of the 
Lagrangian reduces to the boundary part: 



SL 



S(x 3 )d a ( PK a Sz 



K \ 



(4.26) 



The whole dynamics of the system "matter + gravity" 
is, therefore, equivalent to the above equation. Similarly 
as in equation (4.11), we can use the definition of -k^ uk 
and express it in terms of contravariant tensor density 
n^ v obtaining 



' a Ui his — " v^vv ■ (4.27) 

Hence field equations can be written in the following way: 

SL = 8 K (n^SA^) + S(x 3 )d a { PK a Sz K ) . (4.28) 

We complete this Section with the Noether Theorem 
for the energy-momentum tensor (4.18) which, due to 



Bianchi identities (3.9), provides the necessary consis- 
tency condition for the Einstein equations (4.22) or, 
equivalently, (4.24). In fact, due to the regular part 
(4.21) of Einstein equations, Bianchi identity (3.9) re- 
duces to: V a G a b = 0. Hence, the following identity is 
necessary and sufficient for the consistency of the singu- 
lar part (4.22) of Einstein equations: 

Theorem 1. (Noether): For any field configura- 
tion (z K {x)) satisfying the matter dynamics (4.25), the 
energy-momentum tensor (4.18) carried by this configu- 
ration satisfies the following identity: 



V Q r a b = 



(4.29) 



The proof of the Noether identity is given in the next 
Section, just after the Bclinfante-Rosenfeld theorem. 



V. HAMILTONIAN STRUCTURE OF THE 
THEORY 

The above form of field equations is analogous to the 
Lagrangian form of the dynamics in theoretical mechan- 
ics, which may be written as follows: 



SL(q,q) = -jripSq) = pSq + pSq 



(5.1) 



and contains relation between momenta and velocities: 

dL 



P = 



dq 



as well as Newton equations: 
• - 

P dq 

This formula is a starting point of derivation of the 
Hamiltonian form of the dynamics. It is sufficient to per- 
form Lagrange transformation between p and q, putting: 

pSq = 5{pq) - q6p , 

and move the total derivative S(pq) to the left hand side of 
the equation (5.1). This way we obtain the Hamiltonian 
formula: 



-5H(p, q) = pSq — qSp 



(5.2) 



where we have put H(p, q) := pq — L. This formula is 
equivalent to the Hamiltonian form of the equations of 
motion: 



dH 

dp 



P 



dH 



In order to derive the Hamiltonian formulation of the 
field theory, we perform a similar Legendre transforma- 
tion between time derivatives of the fields and corre- 
sponding momenta. For this purpose we have to fix a 
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(3+l)-decomposition of the space-time M. This way the 
theory becomes a Hamiltonian system, with the space of 
Cauchy data on each of the three-dimensional surfaces 
S t := {t = const.} playing role of an infinite-dimensional 
phase space. Unlike in case of the classical mechanics, 
the dynamics of such a system is not uniquely defined, 
unless we control also boundary data for the field in an 
appropriate way. 

In the present paper we consider the case of an asymp- 
totically flat space-time and assume that also the leaves 
S t of our (3+l)-decomposition are asymptotically flat at 
infinity. To keep control over 2-dimensional surface in- 
tegrals at spatial infinity, we first consider dynamics of 
our "matter + gravity" system in a finite world tube T, 
whose boundary carries a non-degenerate metric of sig- 
nature (— , +, +). At the end of the day we shall shift the 
boundary dT of the tube to space-infinity We assume 
that the tube contains the surface S together with our 
matter travelling over it. 

To simplify calculations we choose coordinate system 
adapted to this (3+l)-decomposition. This means that 
the time variable t = x° is constant on three-dimensional 
surfaces of this foliation. We assume that these surfaces 
are space-like. To obtain Hamiltonian formulation of 
our theory we shall simply integrate equation (4.26) (or, 
equivalently, (4.28)) over a finite piece V of the Cauchy 
surface C C M and then perform Lcgcndrc transforma- 
tion between time derivatives and the corresponding mo- 
menta. 

Denoting by V := TflC the portion of the Cauchy 
hypersurface C which is contained in the tube T, we thus 
integrate (4.28) over the finite volume V C C and keep 
surface integrals on the boundary dV of V. They will 
produce the ADM mass as the Hamiltonian of the total 
"matter + gravity" system when we pass to infinity with 
dV = C n dT. Because our approach is geometric and 
does not depend upon the choice of coordinate system, 
we may further simplify our calculations using coordinate 
x 3 adapted to both S and to the boundary dT of the 
tube. We thus assume that x 3 is constant on both these 
surfaces. 

Integrating (4.28) over the volume V we thus obtain: 

6 [ L= f d K {^ V SA; V ) + f S(x 3 )d a ( PK a Sz K ) 
Jv Jv Jv 

(5.3) 

= / (^SA%y+f WA^+f ( PK Hz K ) 
Jv Jdv Jvns 

(5.4) 

where by "dot" we denote the time derivative. In this 
formula we have skipped the two-dimensional divergen- 
cies which vanish when integrated over surfaces dV and 
Vn S (sec fig. 1 below): 




To further simplify our formalism, we denote by px '■= 
Pk° the time- like component of the momentum canoni- 
cally conjugate to the field variable z K . Now we perform 
the Legendre transformation in material variables: 

( PK 5z K ) ' = PK 8z K - z K 5 PK + S { PK z K ) . (5.5) 

The last term, put on the left-hand side of (5.3), meets 
the matter Lagrangian and produces the matter Hamil- 
tonian (with minus sign) , according to formula: 

L mM - p K z K - L nwt - p K °z K o = -T° , (5.6) 

where the canonical energy-momentum tensor is defined 
as follows: 

T\ := pk»z k v - 5? L mat . 

Even if singular, the above quantity satisfies the stan- 
dard Rosenfcld-Bclinfante identity (cf. [11]), which states 
that the canonical energy-momentum tensor is equal 
(modulo a minus sign, due to the convention used here) 
to the symmetric energy- momentum tensor (4.18). 

Theorem 2. (Rosenfeld-Belinfante): Symmetric and 
canonical energy-momentum tensors are, essentially, the 
same. More precisely, the following identity holds: 

T\ = -r^ A g Xv , (5.7) 

or, equivalently, 

T\ = -r ac g cb , (5.8) 

because both the transversal parts: T^„ and t^ v vanish 
identically from the definition. 

Proof. We remember that L mat is a scalar density on S 
and, therefore, may be written as L mat = ^dct g a b A, 
where A is a scalar function, depending exclusively upon 
quantities (z K ; z K a ',9ab)- But the only way to produce a 
scalar from the partial derivatives z K a is to take the fol- 
lowing combination: F KN := z K a z N b g ab . We conclude 
that 

L mat = y/drtg~ b A(z K ,F KN ) . 

This implies identity (5.7) in a straightforward way 
(cf. also [2]). □ 
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Now, we are ready for the proof of the Noether theorem 
(4.29): 

Proof. The invariant character of the matter Lagrangian 
i mat = L mBt (z K ; z K a ; g a b) means that, for any vector field 
X on 5, dragging the arguments (z K ; z K a ;g a b) along X 
produces the same effect on the Lagrangian that dragging 
it directly as a scalar density. Choosing any coordinates 
(x a ) on S and choosing X — d a we obtain, therefore, the 
following identity: 



A. Legendre transformation in the gravitational 
degrees of freedom 

To perform also Legendre transformation in gravita- 
tional degrees of freedom — analogous to transformation 
(5.5) in material degrees of freedom — we follow here a 
method proposed by one of us (see [10]). We show in the 
Appendix that, after the transformation, formula (5.13) 
assumes the following form: 



d a L„ 



dL„ 



0L„ 



8z K 



K 



dz K b 



dz 

+ lr cd d a9cd 



..K 



a+db 



f dL„ 



.,K 



\dz K b 



(5.9) 



where we have used the symmetry of second derivatives: 
d a z K b — dbZ K a and the definition (4.18) of the symmet- 
ric energy-momentum tensor r. Putting now the last 
two terms on the left-hand side and using the Rosenfeld- 
Bclinfantc theorem we obtain: 



d b T b a --r cd d a g cd 



dz K 



dL„ 



dz K h 



Z K a ■ (5.10) 



It may be easily checked that the left hand side is pre- 
cisely the covariant divergence VfcT b a on S (remember 
that r is not a tensor but the tensor density!). Hence, 
we obtain a kinematic identity, fulfilled by arbitrary field 
configurations, not only those fulfilling field equations: 



dz K 

This completes the proof. 



-db 



dz 



K , 



K 



(5.11) 
□ 



Now, integrating equation (5.6) over the two- 
dimensional region V n S, we obtain the "material" part 
of the total ("matter + gravity") Hamiltonian: 

/ L^-p K z K = - [ T% = f r% = / T° , 
Jvns Jvns Jvns Jv 

(5.12) 

where T' 1 " := <5(x 3 )r M!y . Hence, the Legendre transfor- 
mation in material degrees of freedom gives us the fol- 
lowing formula: 

5 f (T°o + L graw )= / (n^SA%y+ f tt^SA^ 

JV JV JdV 

Jv 



+ I (tt k 5z k - z K 6ir K ) 



where we denote 



(5.13) 



(5.14) 



Here, the matter degrees of freedom are already described 
in the Hamiltonian picture (with the matter Hamiltonian 
"— T°o" on the left hand side) and the gravitational de- 
grees of freedom still remain on the Lagrangian level. 



(5.15) 



1 



16-7T J V 

TT K 5Z K 



(P kl 5g kl - g kl SP kl ) +-^J dv (^ Sa - * SX ) 



z K 6w K ) 



lD7T Jgy 



ab 



(5.16) 



Here, we have introduced the following notation: P kl 
denotes the external curvature of S written in the 
ADM form, i.e., given by the equations: 



P kl :=y/dctg mn (Kg kl 
1 

K kl := 



V\9 



001 



^0 

kl 



K kl ) , 

1 



(5.17) 



,001 



1 kl ) 



and g kl stands for three-dimensional contravariant met- 
ric, invariant to the metric g k i induced on the Cauchy 
surface V. Similarly, Q ab denotes the external curvature 2 
of the tube dT written in the ADM form, i.e., three- 
dimensional tensor density given by equations similar to 
(3.1) and (3.2): 



Q ab ■ — V\ dctg cc i| (Lg ab - L ab ) , 

Lab ■ = 



1 ab — 



A 3 

A ab i 



(5.18) 
(5.19) 



and g ab is a three-dimensional contraindicant metric on 
the tube dT, invariant to the induced metric g a b- More- 
over, A = v / deTgAB denotes the two-dimensional volume 
form on dV, whereas 



,30 



a := arsinh 



(5.20) 



is the hyperbolic angle between the Cauchy surface V 
and the tube dT. 

The formula (5.15) has been derived in the paper [3] for 
a wide class of matter Lagrangian (including also gauge 
fields), but only for models with the continuous matter 
distribution. Now we have proved its validity also in case 



2 We use the symbol Q for denoting external curvature of the world 
tube dT to distinguish it from external curvature of the shell S, 
which is denoted by Q. 
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of a singular matter, concentrated on a two-dimensional 
shell S, whose internal metric is non-degenerate and car- 
ries signature (—,+,+). 

Observe that the first term on the left hand side of 
(5.15) vanishes identically due to Einstein equations 3 
which means that the volume part of the total "grav- 
ity + matter" energy vanishes identically. This does not 
mean that the energy vanishes, because there is also an- 
other, surface contribution to the Hamiltonian. Indeed, 
a detailed analysis (see [10]) shows that we are not al- 
lowed to control freely the tube external curvature Q ah 
in the last integral of (5.16) because of constraints which 
occur here. The simplest way to overcome this difficulty 
is to perform another Lcgcndre transformation in the ex- 
pression g ab SQ ab = g oa 5Q°° + 2g 0A SQ 0A + g AB 8Q AB . 
Namely, we write: 



9abSQ ab = S {jjAB Q ) - Q Sg AB , 



(5.21) 



and put the complete derivative 5 (gABQ AB ) on the left 
hand side of (5.15). This way we obtain the "quasi-local" 
(i.e., assigned to the two-surface dV) Hamiltonian of the 
system. Finally, we have: 



-8M 



dV 



+ 



±-Jjx5a-aSx) (5.22) 

/ (n K Sz K - z K 5ir K ) 
Jv 

/ 5oo<5Q 00 + 2g aA 6Q 0A Q AB 5g AB , 

(5.23) 



where the quasi-local Hamiltonian (mass) assigned to the 
two-dimensional surface dV is defined as follows: 



M dv = E (dV) - -L f (Q m g o) 
!6tt Jdv 



(5.24) 



and the additive constant E (dV) is arbitrary. It turns 
out (cf. [10]) that it may be chosen in such a way that the 
quasi-local mass vanishes for the flat Minkowski space 
initial data . The final step of our derivation consists 
in shifting the tube T to infinity. For this purpose we 
limit ourselves to the asymptotically flat case and assume 
that the limiting case of V is equal to an asymptotically 
flat Cauchy three-surface C, i.e., V —* C. It may be 
easily checked that the limit of our quasi-local mass is 



3 The quantity (g° - 87rT° ) in (5.15) is often denoted by NH + 
N k Hk, where H and are the scalar and the vector constraints 
respectively. 

4 Actually, the quasi-local Hamiltonian formula (5.22) needs some 
minor improvements, sec [10], but it is sufficient for global pur- 
poses of the present paper and this is why we do not discuss 
them here. 



equal to the ADM global mass, i.e., that Aigv — ► -Madm- 
Moreover, all the surface terms on the right hand side 
of our Hamiltonian formula (5.22)-(5.23) vanish. This 
way we obtain the following global Hamiltonian formula, 
fully analogous to the mechanical formula (5.2), with the 
ADM mass Madm playing role of the total hamiltonian 
of the "matter + gravity" system: 

SM ADM = ^J (P kl Sg k i - 9kiSP kl ) 

(n K Sz K - z k Stt k ) . (5.25) 

The phase space of our system is thus described by the 
gravitational variables (P ,gui) on the Cauchy three- 
surface C, and by the material variables (wk,z k ) con- 
centrated on the two-surface S C\C, cf. (5.14). 

VI. CONSTRAINTS 

As usual, Gauss-Codazzi equations imply con- 
straints, which must be fulfilled by the Cauchy data 
(P kl , gki, ttk, z k ) on a three-surface C. Outside of the 
shell these are standard, vacuum constraints. In this Sec- 
tion we are going to derive the complete description of 
constraints, valid not only for the regular but also for 
the singular part of the data. We denote by g kl the 
three-dimensional metric inverse to the metric gki and 

(3) 



put 7 := ^/dctgkl■ By R denote the three-dimensional 
scalar curvature of gki, P ■= P kl gu and "|" is the three- 
dimensional covariant derivative with respect to g^i- 

Outside of S Gauss-Codazzi equations relate the com- 
ponents G°f_i of the Einstein tensor density with the 
Cauchy data in the standard way. The spatial part of 
these constraints, tangent to Ct, reads as 



g"i = -Pi 



i \k 



(6.1) 



and the time- like part - normal to Ct, as 

20V = -7 R + (p kl Pki - \P 2 ^J ^ , (6.2) 

where n denotes the future oriented, orthonormal vector 
to the Cauchy surface C t : 

(6.3) 



w = 



9" 



Vacuum Einstein equations outside and inside of S imply 
vanishing of the regular part of Q ^. Hence, the regular 
part of the vector constraint on C reads: 



reg (^ fc | fe ) = , 



(6.4) 



whereas the regular part of the scalar constraint reduces 
to: 

reg (7 R - (p kl Pki - \P^J = . (6.5) 
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In the case of discontinuous matter, the above standard 
constraint shall be completed by their singular part, with 
support on the intersection St — C t n S in the following 
way. Singular part of three dimensional derivatives of 
the ADM momentum P k i consists of derivatives in the 
direction of x 3 : 

sing(P ; fe | fe ) = sing^P, 3 ) = <5(:r 3 )[P 3 ] , (6.6) 

so the full Gauss-Codazzi equations (6.1) take the form 

Q ! = -sing(P ; fe | fc ) = -[P ; 3 ]<*(x 3 ) . (6.7) 

Components of the ADM momentum P kl are regular, 
hence singular part of the term (P kl P k i — 5-P 2 ) vanishes. 
Singular part of the three-dimensional scalar curvature 
consists of derivatives in the direction of x 3 of the (three- 
dimensional) connection coefficients: 

(3) 

sing( J R) = sing(a 3 (r 3 fe ,^-r™,g 3i )) 
= 8{x*)[Tl l ~g kl -T2 l ? 1 ] , 

and expression in the square brackets may be reduced to 
the following term 



7 [rli~9 kl - rzf 1 ] = 



#3 7V<? 33 



75 3fc 



(6.8) 



because derivatives tangent to S are continuous. But ex- 
pression in square brackets is equal to the external cur- 
vature scalar k for the two-dimensional surface St C Ct: 



jk = -d k 



7ff 3fc 



The tangent to St part of Einstein equations gives the 
following: 



G\ = 8irS{x 3 )T° k , 



(6.11) 



which, due to (6.1) and (6.7), implies the following con- 
straints: 



[P 3 k ] = -8nT° h 



(6.12) 



The remaining null tangent part of Einstein equations 
reads: 

where E := T° ^ = T° fl n^d(x 3 ) = 8(x 3 )e, describes 
matter density on C, and s = r M n M , is a surface en- 
ergy density on S. Hence sing(£° M n M ) = 8tte, and scalar 
constraint takes the form 



8ire + [Xk] = 



(6.13) 



(6.9) 



Finally, summing up the regular part of the constraints 
(i.e., equations (6.4) or (6.5) respectively) together with 
their singular parts (i.e., equations (6.6) or (6.10) respec- 
tively) we may finally write down both constraints in its 
distributional form: 

Pl k \ k = -8ttt\ S(x 3 ) , (6.14) 
7 R - [p kl Pki - \p 2 ^ = -167T6 6{x 3 ) , (6.15) 

where the matter momentum r° k and energy s must be 
expressed in terms of material variables (pk, z k ) via the 
matter constitutive equations (4.18). 



VII. CONCLUSIONS 



This implies that 



(3) 



sing Rj = 2 1 ^~g 33 [k]8{x i ) = 2[Xk}S(x 3 ) . 

Finally the total space-like Gauss-Codazzi equation (6.2) 
takes the following form: 



n » 



2g\n 



(3) 

-sing I 7 R 



P H Pkl - 



1 



2 J j J (6.10) 
= -2[Xk]d(x 3 ) . 



Equations (6.7) and (6.10) give a generalization (in the 
sense of distributions) of the usual vacuum constraints 
(vector and scalar respectively). 

Now, we will show how the distributional matter lo- 
cated on St determines the four surface quantities [P 3 k ] 
and [Xk], entering into the singular part of the con- 
straints. The tangent (to S) part of splits into the 
two-dimensional part tangent to S t and the transversal 
part. 



We have proved that the general scheme, used in [10] 
and [3] to describe any continuous, self-gravitating mat- 
ter, may also be extended to the singular matter, concen- 
trated on a 2-dimensional shell. The main result of this 
paper: the Hamiltonian of the complete "matter + grav- 
ity" system is always equal numerically to the ADM mass 
at infinity, similarly as in continuous models. The above 
structure was used in [12] to derive the canonical formu- 
lation of a spherically symmetric dust shell. Recently, 
it was proved that this result may be easily extended far 
beyond the dust case (see [13]). We stress, that the ADM 
mass generates the Hamiltonian evolution of the system 
with respect to the asymptotic time variable at space 
infinity, whereas the local redefinition of the Cauchy sur- 
face C t ■= {t = const.} surface, which does not change 
it at infinity, is merely a gauge transformation. If we 
want to use another time (i. e. the Minkowski time in- 
side the shell or the shell's proper time) there will be 
another Hamiltonian generating the evolution in the new 
parametrization. These issues were thoroughly discussed 
in [14]. 
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Appendix 



Next, we use the following, obvious identity 



To prove formula (5.15) via Legendre transformation 
in gravitational degrees of freedom take metricity condi- 
tions (3.4)-(3.5) (equivalent to equations (2.12)-(2.13)) 
for the connection T on the surface dT = {x 3 = const.} 
and plug them into expression n^dA 3 = 7r 33 <5^4§ 3 + 
2TT 3a SAl a +TT ab SAl b . A straightforward calculation leads 
to the following result: 



n^SA 3 ^ 



Ion 



+d a n 33 6 



-?>a 



(7.1) 



where Q ab denotes the external curvature of the tube dT 
written in the ADM form, i.e., a three-dimensional tensor 
density given by equations similar to (3.1) and (3.2): 



Q ab : = ^\det g cd \(Lg ab 
Li„h ■ = - 



L 



ab \ 



,33 



.,:>>:>> 



A 3 

A ab ) 



(7.2) 
(7.3) 



and g ab is a three-dimensional contraindicant metric on 
the tube dT, inverse to the induced metric g a b- 

Replacing now x 3 by x°, we obtain analogous metricity 
conditions on the surface V = {x° = const.}: 



A 



1 



oo l dfe7r 



A 

A 0k 



A 3 kl n kl ) , (7.4) 

^ (d k n™ + 2A° kl n 01 ) . (7.5) 

Plugging them into expression ir^SA^ = tt 00 SAq + 
2-7r 0fe <L4 3 lfe + TT kl SA° kl we obtain immediately the follow- 
ing identity: 



T ^SA° 



1 

16tt' 



J)k 



Zoo > (7-6) 



where P kl denotes the external curvature of £ written in 
the ADM form, i.e., given by the equations: 

P kl :=Vdct<7 mn (Kg kL -K k 
J^ki ■— - 



(7.7) 



V\9 



oo 



:i m — - 



V\9 



oo 



4° 

■ A kl ) 



and g kl stands for three-dimensional contravariant met- 
ric, inverse to the metric g k i induced on the Cauchy sur- 
face V. 

Using these results and skipping the two-dimensional 
divergencies which vanish after integration over dV, we 
may rewrite the gravitational part of (5.4) in the follow- 
ing way: 

j y (^SA%y + n^SA^ = -^J v {9kiSP kl ) 

~j 9ab6Q ab 
lOTT Jgy 



y/\TT°°TT 33 \ 

(7.9) 



Denote 



We have 



r 30 



„30 



(7.10) 



2Vk 00 7T 33 | = 



A 



where A := \J det gAB ■ Hence, we obtain: 



7T 5 



7T° 3 \ „,/^ 3 °\ 1 



-00 



+ 7T 33 5 



IT"" J \TT' 

where a := arcsinh(q) and, consequently 



(7.11) 
(7.12) 



/ v (^)-+/ w ^ = -isr/>^)- 

/ 9ab5Q ab + ±- I (XSaY . 

lOTT Jgy 8tT Jgy 

Now we perform Legendre transformation between time 
derivatives and the corresponding canonical momenta. 
This transformation is preformed both in volume: 



(9kiSP kl )' = (g k iSP kl -P kl Sg kl 



S(g k iP kl 



dV 



„33 < 



(7i 



and on the boundary: (XSa)' = (XSa— aSX)+S(Xa). The 
sum of the two total derivatives which arise here may be 
easily calculated using the same arguments as in paper 
[10]: 

(7.13) 

Moving the first (volume) quantity to the left hand side 
of formula (5.13) and collecting it with the gravitational 
part of the Lagrangian, we obtain: 

(7.14) 

which may be treated as the "volume part of the gravita- 
tional Hamiltonian" . It meets the "matter Hamiltonian" 
(5.12). Their sum (the "volume part of the total Hamil- 
tonian" ) vanishes identically as a consequence of Einstein 
equations. This completes the proof of formula (5.15). 
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